Q- Value and Half-Lives for the Double-Beta-Decay Nuclide ^^°Pd 
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The ^^"Pd double-beta decay Q-value was measured with the Penning-trap mass spectrometer 
ISOLTRAP to be Q = 2017.85 (64) keV. This value shifted by 14keV compared to the literature 
value and is 17 times more precise, resulting in new phase-space factors for the two-neutrino and 
neutrinoless decay modes. In addition a new set of the relevant matrix elements has been calculated. 
The expected half-life of the two-neutrino mode was reevaluated as 1.5 (6) x 10^" yr. With its high 
natural abundance, the new results reveal ^'^"Pd to be an excellent candidate for double-beta decay 
studies. 

PACS numbers: 07.75.+h, 14.60.Pq, 23.40.Bw, 32.10.Bi 



The recent results on neutrino oscillations have 
revolutionized our understanding of the role played by 
neutrinos in particle physics and cosmology, in particular 
by proving that neutrinos have a finite mass. In the quest 
for a detailed understanding of the neutrino itself, the 
rare process of double-beta decay offers the most promis- 
ing opportunity to probe the neutrino character and to 
constrain the neutrino mass 0, 0| • In contrast to neutrino 
oscillations, which violate the individual flavor-lepton 
number while conserving the total lepton number, the 
process of neutrinoless double-beta decay (Oz^/3/3-decay) 
violates total lepton number and is as such, forbidden by 
the Standard Model of particle physics. Moreover, un- 
like the observed neutrino-accompanied double-beta de- 
cay (2z//3/3-decay) process, the Oz//3/3-decay process would 
imply that the neutrino is a Majorana particle, i.e., its 
own antiparticle. 

While the decay spectrum of the 2i^/3/3-decay is contin- 
uous, the experimental signal of the Oz//3/3-decay repre- 
sents the sum energy of the two electrons at the decay 
Q-value. The expected half-life of the Oi//3/?-decay is ex- 
tremely long and hence, very small event rates are ex- 
pected. In addition, a high accuracy (below 1 keV) is 
desirable to properly identify the signal with respect to 
background. Furthermore, a well-known Q-value allows 
a precise determination of the phase space of the half- 
lives of the double-beta decay modes. 
The decay rates of both decay modes are strong func- 
tions of the Q-value. Oi^/3/3-decay scales with the fifth 



power of the Q-value and 2i^/3/3-decay with the eleventh 
power. Eleven nuchdes (^^Ca, '^''Ge, ^^Se, ^^Zr, ^°°Mo, 
"Opd, ii^cd, i24sn, i30Te, "exe, and i^ONd 0) are con- 
sidered as potential Oz//3/3-decay candidates, having high 
enough Q-values (above 2 MeV) for relatively short half- 
lives, which are accessible within the experimental limits 
and natural abundances for feasible experiments. Among 
these, ^^'^Pd has the second highest abundance, making 
it a particularly interesting case. Furthermore, ^^°Pd 
is an excellent test candidate for the "single-state dom- 
inance hypothesis" (SSDH) [§-[1^, which predicts that 
the double-beta decay rate is dominated by a virtual two- 
step transition through just one single intermediate 1"*" 
state. Nevertheless, the experimental knowledge about 
the double-beta decay of ^^"Pd is still very limited. Re- 
cently, first experimental limits on excited state transi- 
tions were published However, ^^'^Pd has the highest 
Q-value uncertainty of all candidates mentioned above. 
In this Letter, we report the first direct mass comparison 
between ^^°Pd and the double-beta decay daughter nu- 
clide ^^"Cd. The new Q-value is shifted by almost 14keV 
to the literature value [12| and has a 17 times smaller un- 
certainty. Furthermore, the absolute masses of ^^°Pd and 
^^"^Cd were measured with high precision within the same 
measurement campaign. We also report new calculations 
for the phase-space factors of the /3~/3~-double-beta de- 
cay modes based on the new Q-value. In addition, the rel- 
evant nuclear matrix elements have been calculted as well 
as the resulting half-lives of the neutrino-accompanied 
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FIG. 1. Scheme of ISOLTRAP as used for the present mea- 
surements. The atoms are ionized with the ablation ion 
source, trapped in the first trap for purification and trans- 
ferred to the second trap to measure their cyclotron frequency. 



FIG. 2. Time-of-flight resonance as a function of the excita- 
tion frequency for ^^^Pd"*" for the Ramsey excitation scheme 
50-600-50 ms with a fit of the expected theoretical line shape 
(solid line) [3. 



and neutrinoless double-beta decay. 
The measurements were performed with the Penning- 
trap mass spectrometer ISOLTRAP [13.] at the CERN- 
ISOLDE facility. ISOLTRAP determines the masses by 
measuring the cyclotron frequency Vc = qB / {2iTm) of the 
corresponding ions with mass m and charge <? in a mag- 
netic field B. For the present off-line-experiment, a laser 
ablation ion source for the production of carbon-cluster 
reference-mass ions IJ] was modified to deliver singly 
charged ^^°Pd and -'^^'^Cd ions. In short, a frequency- 
doubled Nd:YAG laser was focused with a typical pulse 
power density of 5 • lO'^ W/cm^ on either a carbon, a pal- 
ladium or a cadmium foil for desorption and ionization. 
All three targets were mounted on a rotary sample holder 
for rapid element selection. 

After the ablation from the sample, the ions were guided 
to the preparation Penning trap (see FigH]), where they 
were centered and purified by buffer-gas cooling [l5j. 
They were then transferred to the precision Penning trap 
where the time-of-flight ion-cyclotron resonance (TOF- 
ICR) detection technique [16| was used to determine their 
true cyclotron frequency Vc- Two different measurement 
sets using two different Ramsey excitation schemes [TtI - 
[igj were performed to increase the accuracy. The first ex- 
citation scheme (30-840-30 ms) consisted of two rf-pulses 
of 30 ms separated by a 840 ms waiting period. For each 
resonance curve (Fig.[2|), data were taken for about 20 
minutes. In total, 17 such resonances of ^^°Cd+ and 
16 resonances of ^^'^Pd"'" were recorded alternately. For 
the second set of measurements, a 50-600-50 ms Ram- 
sey excitation scheme was applied with the same dura- 
tion of data-taking per resonance and the same num- 
ber of resonances. The Q- value of the double-beta decay 
of the mother nuclide with mass nim to the daughter 
nuclide with mass is given by the mass difference 



Q = m„i — rud, which in turn can be written as a func- 
tion of the frequency ratio r = Vd/vm and the electron 
mass me'- 

Q = m„i -nid^ 1 — ^ 1 ) {"^d - me). (1) 

V t^m } 

Figure[3] shows the resonance frequencies for the 50-600- 
50 ms Ramsey excitation scheme as a function of time of 
the measurement campaign with typical statistical uncer- 
tainties on the order of cr{vc)/vc = 1 • 10^*. In addition, 
different systematic uncertainties have to be taken into 
account: the presence of contaminating ions, the time- 
dependent magnetic-field changes, the mass-dependent 
systematic effect, and the residual systematic uncertainty 
of ISOLTRAP [20| . A countrate-class analysis applied to 
the individual frequency measurements confirmed that 
there were no contaminating ions or space-charge effects 
present The "°Cd+ and "Opd+ data from Fig.[3] 

follow the same trend, mainly given by a variation of 
the magnetic field strength in time. Following Bradley 
et al. 21|, the frequency ratio r and thus the mass ratio 
between the two ion species is deduced by fitting simul- 
taneously a pair of polynomials 



Mt) ^ Cq + Ci ■ t + C2 ■ + ... + Cn ■ t" 



(2) 



Mt) = r ■ {cq + ci ■ t + C2 ■ + ... + Cn ■ i") (3) 

to the datasets of "°Pd+ and "°Cd+. Here, the fre- 
quency ratio r entered as one of the fit parameters, while 
the other fit parameters represent the behavior of the 
magnetic field strength in time. Polynomials (Eq. ([2]) and 
Eq. ([3])) of different orders were tested with smallest 
found for n = 3. The procedure was performed sepa- 
rately for both sets of Ramsey schemes. In Fig.O the 
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FIG. 3. Cyclotron frequencies of ^^"Cd^ (circles, top) and of 
iiOp^+ j^gqQares, bottom) using the Ramsey excitation pat- 
tern 50-600-50 ms. The solid lines are simultaneous fits to 
both isobars (see text). 



result of the 50-600-50 ms data is shown as an example, 
resulting in a reduced of 1.08. In case of the 50-600- 
50 ms data, the of the polynomial was 1.10. The mass- 
dependent effect is more than one order of magnitude 
smaller than the achieved statistical uncertainty. Finally, 
an additional relative systematic uncertainty is taken into 
account for the slight differences between the production 
and ion-guiding parameters of ^-^^Cd"*" data and -^^°Pd+. 
Although this uncertainty is expected to lie well below 
the statistical uncertainty for mass doublets, the rela- 
tive residual systematic uncertainty of ISOLTRAP for 
frequency ratios of (5r/r — 8 • 10~^ was taken as an up- 
per limit Table U lists the obtained frequency ratios 
and Q-values as derived from Eq. ([ij of both measure- 
ments, their weighted averages, and the Q- value based on 
AME2003 [ll|. The resuhs from the two data sets agree 
with each other within 0.55 a and the uncertainty of the 
weighted average of both measurements is 0.64 keV. This 
is a reduction by more than a factor of 17 from the liter- 
ature value based on the AME2003 [12j. 
Within the same measurement campaign, not only the 
Q-value, but also the absolute masses of ^^°Cd and 
^^°Pd have been determined using carbon-clusters as 
reference ions. Their uncertainties were reduced by a 
factor of 1.7 and 5, respectively. While the mass ex- 
cess ME = -90348.84 (1.30) keV of "°Cd is in good 
agreement with the literature value, the mass excess 
ME = -88333.23 (1.34) keV of "°Pd deviates by more 
than 20 keV JJj . Due to the systematic mass-dependent 
effect, which does not cancel using reference ions of dif- 
ferent masses, the extracted Q-value of 2015.6(1.9) has 
a larger uncertainty, but agrees well with the direct mea- 
surement of the Q- value. 

The phase space of the 2i//3(3-decsLy and of the 
Oz^/3/3-decay and thus their half-lives are affected by the 



TABLE I. Results of the evaluation of the Q-value measure- 
ments of the ""Pd and ""Cd mass doublet. Second column: 
cyclotron frequency ratios; third column: Q-values with un- 
certainties, where the residual systematic uncertainty domi- 
nates the value by more than 70%. First and second line: 
individual Ramsey data sets; third line: weighted avera ge o f 
both data sets; fourth line: Q- value based on AME2003 [1^ . 



Data 


Ratio r 


Q-valuc / keV 


30-840-30 ms 


1.0000197131(89) 


2018.09(90) 


50-600-50 ms 


1.0000197081(89) 


2017.60(90) 


Weighted Average 


1.0000197106(63) 


2017.85(64) 


AME2003 




2004(11) 



new measurement of the Q-value. Half-lives and phase- 
space factors G are linked via 

{T^J^y ^G^''iQpp,Z)g\\m,Mlr^\' and (4) 
(TO;,)-i = G''^{Q,,,Z)g\ I M'-^ f\ ^ f , (5) 

Tile 

where mg is the mass of the electron, g^^ is the axial- 
vector coupling constant, and (m^) is the effective Ma- 
jorana neutrino mass. M^^rp and Af°'' are the matrix 
elements of the corresponding decay modes. 
The new Q- value is almost 14keV higher and thus, in- 
creases the phase-space of the transitions and shortens 
the expected half-lives of both decay modes. A general 
formulation of the calculation of the phase-space factors 
was given by Doi et al. [2^ . However, in contrast to pre- 
vious calculations, where an approximate expression for 
the electron wave function at the nucleus was used, the 
phase-space factors reported here are calculated by a new 
method [23*1 using exact Dirac electron wave functions 
and including electron screening by the electron cloud. 
The wave functions are obtained by numerically solving 
2J] the Dirac equation with a potential for which the 



electron screening is taken into account through the nu- 
merically obtained Thomas-Fermi function [25|. A nu- 
clear radius of i? = tqA^^^ with rg = 1.2 fm was assumed 
for these calculations. The possible sources of uncertain- 
ties for both decay modes are the Q-value, the electron 
screening, the nuclear radius, and in case of the neutrino- 
accompanied mode, the closure energy. The uncertainties 
coming from the Q-value, which is dominating the total 
uncertainty in case of the old Q-value, could be reduced 
by more than an order of magnitude due to the more 
precise measurements. A comparison between the space 
factors based on the new Q-value and the Q-value from 
AME2003 and their uncertainties is shown in tablelTTl 
Using the single-state-dominance hypothesis (SSDH) 
ilOj , which works well for the nearby nuclei lOOMo and 
116Cd, the SSDH-half-life was calculated to be T^^^ = 
1.5 (6) X 10^'^yr, which is in good agreement with pre- 
viously predicted values [1, M, [23 ■ The uncertainty of 
the half-life was evaluated considering the uncertainties 
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TABLE II. Comparison of the phase-space factors for 
2!//3/3-decay and 0!^/3/3-decay based on Q-values determined 
from the AME2003 mass values [12| and the new ISOLTRAP 
values. The values in the first brackets are the uncertainties 
caused by the uncertainty of the Q-value, while the values in 
the outer right brackets are the total uncertainties. 



Q,3/3 G^" / yr~' 


/ yr-^ 


2004 1.313 (64) (78) X 10""'' 


4.716(23) (33) x 10"'" 


2017.8 1.391 (4) (19) x 10"''=' 


4.824(2) (12) X IQ-''' 



of the /3-decay half-life and the minimal and maximal 
value of the /3-decay and electron-capture amplitudes of 
the 1+ intermediate state ^^''Ag. In general, the SSDH 
half-life is independent of the value of the axial coupling 
constant 0j since it is included in the log /i- values. 
Nevertheless, an axial coupling constant of a free nucleon 
{gA = 1.269) was assumed, in order to calculate the ma- 
trix element for the neutrino-accompanied double-beta 
decay to be M^^j, = 0.263 MeV'^ 

To specify the required sensitivity of future Oi^/3/3-decay 
searches, the nuclear matrix element M^" for ^^'^Pd was 
calculated using two different models: the Quasiparti- 
cle Random Phase Approximation theory (QRPA) [28] 
and the Microscopic Interacting Boson Model (IBM-2) 
[ioj , using the CD-Bonn and Argonne potentials for the 
short-range correlations as in i28j]. The resulting ma- 
trix element using QRPA with the CD-Bonn potential is 
M^" = 6.5 (9) and M^" = 5.9 (8) with the Argonne po- 
tential. With the IBM-2 model one gets M^" = 4.3 (1.3) 
and M"'' = 4.1 (1.2), respectively. We assumed for 
these calculations gA = 1.269 and a nuclear radius of 
R = roA^^^ with tq = 1.2 fm. These matrix elements 
have relatively large values compared to most other iso- 
topes [Uli^. 

A reliable estimation of the half-life of the neutrinoless 
double-beta decay is more complicated than in case of 
the neutrino-accompanied mode. Here, the uncertainty is 
mainly due to the unknown neutrino mass. Using Eq. ^ 
and the weighted average M'^'^ — 5.5 of the four different 
matrix elements, obtained with the QRPA- and IBM- 
2-model and the CD-Bonn and Argonne potentials, the 
half-life can be estimated by 



^1/2 — 



6.8 X 10^'^eV^yr 



(6) 



Assuming the effective neutrino mass in the range from 
1 eV down to 10~'^eV leads to a 0:^/3/3-decay half-life 
between 6.8 x 10^'^ yr and 6.8 x 10^^ yr, respectively. 
Currently, a possible evidence for Oz^/3/3-decay in ^^Ge 
is discussed resulting in an effective Majorana neutrino 
mass of {m„) = 0.32 (3) eV 3l|. To explore this mass 
range, a Oi^/3/3-decay experiment on ^^'^Pd must be sen- 
sitive to a half-life range from — 5 x 10^* yr to 
T^J^ = 1 X 10^^ yr. 



In conclusion, the Q-value of the double-beta transition 
from ^^°Pd to ^^°Cd was measured for the first time by 
high-precision Penning-trap mass spectrometry and re- 
suhed in a value of Q = 2017.85(64) keV. Thus, the un- 
certainty was reduced by a factor of 17 with respect to 
the AME2003 value, revealing a 14keV higher Q-value 
compared to the literature value. In addition, masses 
of ^^°Pd and ^^°Cd were determined, reducing their un- 
certainties significantly. Based on the more accurate Q- 
value, the phase-space factors of the 2i//3 /3-decay and the 
Oj^/3/3-decay were recalculated by a new, more precise 
method. The SSDH half-life calculation leads to an ex- 
pected half-life of the neutrino-accompanied double-beta 
decay of T^^^ = 1.5(6) x lO'^^yr. In case of the neutri- 
noless double-beta decay, the matrix elements have been 
calculated with the help of the QRPA and IBM-2 models, 
resulting in large values compared to most other possi- 
ble double-beta decay nuclides. Thus, in combination 
with its high natural abundance, ^^''Pd becomes a very 
promising candidate for double-beta decay studies and 
for the search for the neutrino mass. 
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